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Abstract 

We present an extension of the Randall-Sundrum model in which, due to sponta- 
neous Lorentz symmetry breaking, graviton mixes with bulk vector fields and becomes 
quasilocalized. The masses of KK modes comprising the four-dimensional graviton 
are naturally exponentially small. This allows to push the Lorentz breaking scale to 
as high as a few tenth of the Planck mass. The model does not contain ghosts or 
tachyons and does not exhibit the van Dam-Veltman-Zakharov discontinuity. The 
gravitational attraction between static point masses becomes gradually weaker with 
increasing of separation and gets replaced by repulsion (antigravity) at exponentially 
large distances. 

1 Introduction 

The discovery of cosmological acceleration [1] indicates the necessity for modification of the 
standard cosmological paradigm. One way to explain the accelerated expansion of the Uni- 
verse is to attribute it to the effect of some unconventional matter such as the cosmological 
constant or quintessence. In this approach one encounters, however, the difficulty of explain- 
ing the small value of the associated energy scale in comparison with other fundamental 
scales, such as the Planck mass or electroweak symmetry breaking scale. It is reasonable 
to consider an alternative possibility that the cosmological constant is exactly zero, but its 
effect is mimicked by modifications of gravitational laws at distance scale of the order of the 
present cosmological horizon. 

A natural framework for modifying gravity at ultra-large distances is provided by models 
with infinite extra spatial dimensions. A number of brane- world models have been proposed 
[2, 3] where the four- dimensional graviton is a superposition of bulk KK modes, and thus, is 
not localized on the brane. Rather, it is quasilocalized, having a non-zero width T (and, pos- 
sibly, mass m) with respect to decay into extra dimensions. In these models gravity is four- 
dimensional at distances smaller than r c ~ min{m _1 ,r -1 } and becomes multi-dimensional 
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at larger scales. The attractive feature of the models with extra-dimensions is that the tiny 
values of m, F are generated without strong fine-tuning of parameters. 

However, Lorentz invariant extra dimensional models with infrared modifications of grav- 
ity face the problems [4, 5, 6] akin to those of four-dimensional massive gravity. Either these 
theories contain ghosts, fields with the wrong sign of kinetic term, or the propagator of gravi- 
ton exhibits the van Dam-Veltman-Zakharov (vDVZ) discontinuity [7] originating from a 
scalar degree of freedom that does not decouple in the massless limit. At the classical level 
the discontinuity might be cured by nonlinear effects [8, 9], but at the quantum level the 
presence of the additional scalar leads to strong coupling [10, 5, 6] at the energy scale at most 
of order [m 2 Mpi) 1 ^ 3 . If m is comparable to the present Hubble parameter of the Universe, 
the theory is in strong coupling regime, and thus looses its predictivity, at distances smaller 
than 1000 km. 

There are indications that these difficulties might be resolved in the Dvali-Gabadadze- 
Porrati (DGP) [3] model. In particular, it was argued that the scale of strong coupling can 
be pushed to higher energies in nontrivial classical backgrounds, provided the structure of 
counterterms in the model is of a special form [9]. However, it is still unclear whether the 
counterterms in the DGP model actually meet the requirements of Ref. [9]. Probably, the 
most disappointing fact about the DGP model is that the branch of the cosmological solutions 
in this theory, which exhibits cosmic acceleration without introduction of the cosmological 
constant, is plagued by ghost-like instabilities [5, 9]. 

In this paper we propose a brane-world model with infrared modification of gravity and 
violation of 4-dimensional Lorentz-invariance. In our model graviton becomes quasilocalized 
due to mixing with vector fields which freely propagate in the bulk. This kind of mixing 
becomes possible once the Lorentz symmetry is spontaneously broken by condensates of the 
vector fields. 

Our approach builds up on the ideas put forward recently in the four-dimensional frame- 
work. It was suggested that the problems of massive gravity could be resolved in models 
incorporating violation of Lorentz invariance [11, 12, 13]. Models with explicit Lorentz vio- 
lation in the gravitational sector were considered in Refs. [11, 12, 13, 14] and it was shown 
that they are equivalent to a class of models with gravity coupled to scalar fields with un- 
usual kinetic terms. In the latter formulation the Lorentz symmetry is broken spontaneously 
by coordinate dependent scalar condensates. Because of nonlinear structure of the scalar 
kinetic term in this kind of models, they possess a cutoff scale A cuto f f = (mMpi) 1 ^ 2 , and can 
be considered only as low energy effective a-models. If the graviton mass m is of the order 
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of the present Hubble parameter, then A cuto ff ~ (O.Olmm)^ 1 , which is phenomenologically 
acceptable, but still unnaturally small in comparison with other fundamental scales. It is 
unclear whether these a-models can be extended to complete models making sense above 
the scale A cutoff . 

As a scenario for such a complete model one could envisage a gravitational analog of the 
Higgs mechanism triggered by spontaneous Lorentz symmetry breaking. In Refs. [15, 16, 17, 
18, 19] it was proposed to use coordinate independent vector (in general, tensor) condensates 
for this purpose. However, in four-dimensional theories, graviton does not acquire mass 
within this approach: no gap in the dispersion relation of graviton appears [20, 21]. Moreover, 
the effect of vector condensates on cosmology amounts to nothing but the change of the 
gravitational constant in the cosmological equations [19, 21], and no cosmological acceleration 
is produced. 

The brane-world approach opens up a new way to circumvent the problems encountered 
in four dimensions. In this paper we make the first step along this line: in Sec. 2 we present 
a model in which graviton becomes quasilocalized due to spontaneous Lorentz symmetry 
breaking induced by non-zero VEVs of bulk vector fields, and analyze in the subsequent 
sections the spectrum of linear perturbations about four-dimensionally flat background. In 
Sec. 3 we analyze the tensor perturbations and find that the characteristic mass m c of KK 
modes comprising four-dimensional graviton is naturally exponentially small, which allows 
to take the Lorentz-breaking scale as high as a few tenth of the Planck mass. In Sec. 4 
and Sec. 5 we study vector and scalar perturbations and show that the model contains 
neither tachyons, nor ghosts, and does not exhibit the vDVZ discontinuity. We calculate 
long distance modification of the Newton law and find quite unexpectedly that attraction 
of point masses changes to repulsion at distances of order l/m c . Thus, the model provides 
a realization of antigravity at ultra-large distances. This looks promising in view of the 
possibility to obtain cosmological acceleration at late times. On the other hand, modification 
of gravitational field of static sources in the model faces severe phenomenological constraints. 
We discuss possible ways to get around these constraints in the concluding Sec. 6. 

A number of important issues are left beyond the scope of the present paper. These are, 
e.g., cosmology of the proposed model and structure of quantum corrections. We leave them 
for future investigation. 



3 



2 The model 



We consider a five-dimensional setup with positive tension brane and negative cosmological 
constant in the bulk — the Randall-Sundrum (RS) model [22] . We add to this model three 
bulk vector fields A a M) a = 1,2,3, with quartic potential localized on the brane. The action 
is taken in the following form, 

+ J d'x^g {-a - ^ {r v A^Al + v 2 5 ab ) 2 - ^ [gT A%A h v - V^^) ^ , 

(1) 

where g^ v is the induced metric on the brane, and 

Fmn = &mA n — 8nA m . 

The capital Latin indices M, N take values 0, 1, 2, 3, 5, while the Greek indices /z, z/, . . . run 
from to 3. The action (1) is invariant under global 5*0(3) symmetry with vector fields A a ^ 
belonging to fundamental representation. It is straightforward to check that the potential 
in the second line of Eq. (1) is generic quartic potential invariant under this group. Note 
that the second term in the potential depends only on the traceless part of the matrix 
A a A h v . The parameters v 2 , xj, h\ are assumed to be positive. In addition to the global 
SO (3) symmetry, the bulk part of the action is invariant under U(l) x U(l) x U(l) gauge 
transformations, 

A a M ^A a M + d M a a . 

This Abelian gauge symmetry is broken explicitly at the brane. 

The model has a static solution with AdS metric in the bulk and constant values of the 
vector fields, 

d s 2 = -dz 2 + e-^rj^dx^dx" , (2a) 
At = A a =0, A* = v5t, i = 1,2,3, (2b) 

where k = y^—AnG 5 A/3, and the usual fine-tuning a = —A/k is assumed; the signature 
of the Minkowski metric rj^ is (+,—,—,—). The VEVs (2b) of the vector fields break 
5*0(3) xLorentz symmetry down to 50(3) of spatial rotations accompanied by simultaneous 
rotations in the internal space. The similar pattern of Lorentz symmetry breaking was 
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considered in the four- dimensional context in Refs. [23, 21]. As we will see, this spontaneous 
symmetry breaking provides mixing between the graviton zero mode, present in the pure 
RS case, and the continuum spectrum of KK modes of vector fields. This mixing results in 
quasilocalization of the graviton. 

Let us study the linearized perturbations above the background (2). Imposing the gauge 
conditions = 0, = —1 on the metric, one writes the following decomposition, 

ds 2 = ~dz 2 + (e~ 2fc|z| r^ + h^(x, z))dx fl dx v , 

The analysis is simplified by the fact that the energy- momentum tensor of the vector fields 
in the bulk vanishes to the linear order in perturbations, allowing for fixing the transverse 
traceless gauge in the bulk [22, 24], 

K = = o • (3) 

Here and throughout the paper indices //, u, . . . are raised (lowered) using the metric rf v 
{j]jxv)i and comma denotes derivative. In this gauge the bulk equations for the metric take a 
fairly simple form, 

where u(z) = e~ fc ' 2 ', and prime denotes derivative with respect to z. As to the vector fields, 
the bulk equations read 

-atr + 2ka" + \a*% = 0. (5) 

In deriving this equation we made use of the [^(l)] 3 gauge invariance of the bulk action to 
impose the conditions 

«s = 0, a a % = (6) 

on the vector fields in the bulk. Imposing the two conditions (6) simultaneously is possible, 
as they are compatible on-shell. The coordinate frame where the conditions (3) and (6) are 
satisfied will be referred to below as the bulk frame. 

The boundary conditions on the brane are most easily formulated in the Gauss normal 
(GN) reference frame. In this frame the brane is fixed at z = (the quantities with the 
bar refer to the GN frame), and we assume Z2 symmetry across the brane. Then, for 
general energy-momentum tensor on the brane, the boundary conditions for the metrics 
at z = +0 read [24], 

% v + 2k V = 8ttG 5 (t^ - ^uT^j . (7) 
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The linearized energy-momentum tensor of the vector fields comes from the Lorentz violating 
term in the quadratic action arising from the brane potential in Eq. (1). To the quadratic 
order, the latter term has the form, 

1 



7 (2) 



= - Id x — ^— I h ab 



+ 



x 2 v 4 



h 



III) 



hrr + ~a\ 



2 ~i 



(8) 



From this expression one obtains, 

rpV _ rpV _ n 



2x 2 V f 
—=—6, 



ab 



7 ^ _fc 

hkk H — a-k 



(9a) 
(9b) 



where we introduced x 2 = k\ + x|. Note that the energy- momentum tensor (9) of the vector 
fields violates the weak energy condition. We will show that this property does not lead to 
instabilities. At the same time it is crucial for antigravity emerging at ultra-large distances. 

For the vector fields, the Z 2 symmetry and continuity demand that the 5th component 
vanishes on the brane, a^(z = 0) = 0. The sources for the other components can be read off 
from (8), resulting in the following boundary conditions at z — +0, 

ag' = , (10a) 



-2„.3 



I h ai + -« + <)] -y 



2„,3 



hkk H — Ofc 

v 



(10b) 



The metrics in the bulk frame and GN frame are related by a gauge transformation. The 
form of the latter is restricted by the condition h 55 = h 5/Ji = 0. One obtains 



h^u — h^ v + u (E/j,,u 4~ &u,[m) 4~ ^iV-v 2/uM Tj^ u s , 



;n) 



where the functions e, depend only on x. This transformation corresponds to the following 
change of coordinates, 

1 



z + e 



x/' 



2ku 



Note that in the bulk frame the brane is displaced from the origin along the z-direction, its 
position being described by the brane bending e. The relation between the components of 
the vector fields in the two frames has the form 

a% = a%- ^e, a + d 5 a a , 

-a a V a a 
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where we allowed for the possibility of a [?7(1)] 3 gauge transformation parametrized by the 
functions a a (z,x). Imposing the gauge a§ = 0, which is consistent with the boundary 
conditions on the brane, we obtain, 

a£ = < + v£ l + . ( 12 ) 

where the functions f3 a depend only on x. The set of equations (7), (9)— (11) and (12) results 
in the following boundary conditions, 

A / 2 1 2 \ 

h' 00 + 2kh 00 + 2e fi0 h aa + -< + -e >aa + 6ke + -(3 a >a , (13a) 

+ 2kh 0i + 2e fii = , (13b) 

h'^ + 2khij + 2e,ij = a{ % + -(aj + a{) + ^ + -(/% + 

If 2 1 2 \ 1 

- gtfo- Ufcfc + -aj: + -e ifcfc + -/3*J \ , (13c) 

ag' = , (13d) 
af = xV{/w + I(o? + <) + i e>ai + + /£)} + 2x?t; 3 fo5 aie 

(21 2 \ 

h kk + -Ofe + T e,kk + -P h k ; (13e) 
v k v J 



x 2 y 
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here we introduced the notations 

Ai = 167rG 5 xy , A = 167rG 5 xV . 
In what follows it is convenient to work in the 4-dimensional Fourier representation, 

V, < oc e~ ipoxo+ipx , 

and decompose the perturbations into scalar, vector and tensor modes with respect to ro- 
tations around the three-momentum p. To this end one introduces a three-dimensional 
orthogonal basis consisting of p and two unit vectors, e^ a \ a = 1,2. The latter are used to 
construct a pair of transverse traceless symmetric tensors, 

~ \/2 j j ' j ~ y/2 j j ' 

and a transverse antisymmetric tensor, 

3 \/2 ' 



The tensorial decomposition reads, 

hoo =0i , (14a) 

h i =Vi<\>2 + efVia , (14b) 

hij = (Sij - ^f^j 03 + PiPj(pi + (Vi^ + Pj^) + 4fxia , (14c) 

a a =p a fo + ei a ^ 3a , (14d) 

< = - ^r) + PiPafa + + Pae^) 4a + (p l e[ a) - p a e [ t ] ) 

+ d^ X 2a + M, (14e) 

I3 a = Pa <Ps + ei a Ve a , (14f) 



with p = (pi,P2,P3)- There are two fields x i n the symmetric tensor sector; one field £ in 
the sector of antisymmetric tensors; six fields ip in the vector sector, and eight fields in 
the scalar sector. Note that all fields in (14), except for 8 and i/jq, are functions of the fifth 
coordinate z. 



3 Tensor modes 

Let us study the spectrum of tensor perturbations. To begin with, one notices that the 
antisymmetric field £ completely decouples. It satisfies Eq. (5) with the free boundary 
condition £'(z = 0) = 0, and forms a continuum spectrum of completely delocalized modes 
similar to the modes of a free vector field in RS background [25]. We do not consider these 
modes below. 

For the symmetric tensors one obtains from Eqs. (4), (5) the following equations in the 
bulk, 

l -Xi - 2A; 2 Xl + |^xi = , (15a) 

2 

- x >> + 2kx' 2 -- 2 X2 = U , z>0, (15b) 

where p 2 = pi — p 2 , and polarization index a is omitted to simplify notations. The boundary 
conditions are read off from Eqs. (13c), (13e): 

x [ + 2kxi = \[xi + ^X2^j , (16a) 
X ' 2 = xV ( X i + jW) , z = +0. (16b) 
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The system (15), (16) can be rewritten in the form of an eigenvalue problem, 



-u 2 ( X " - ±k 2 X i) - [(4fc - 2A) X i - 2^x2} 6(z) =m 2 Xl , (17a) 
-u 2 ( X ' 2 ' - 2k X ' 2 ) + [2fi X i + 2u X 2]S(z) =m 2 X 2 , (17b) 



where we denoted m 2 = p 2 , v = 2k 2 v 2 , /x = \f~Xv and introduced X2 = V^2~kUIx2 to make 
the operator Hermitian with the scalar product 




First, let us make sure that there are no tachyonic modes. The operator (17) is the sum 
of the diagonal operator appearing in the pure RS case and a positive semi-definite operator 



Thus, the absence of tachyons in our case is ensured by their absence in the RS model. 

Next, there is no normalizable zero mode. Indeed, if m 2 = 0, then X i = A\e~ 2k][Z \ 
X2 = ^-2 = const (other solutions to Eqs. (15) are exponentially growing). From Eqs. (16) 
one obtains that 2XAi + 2[iA 2 = 0. Thus, A 2 7^ and the zero mode is not normalizable 
with the scalar product (18). 

We now turn to the modes with m 2 > 0. These modes belong to continuum spectrum. 
Introducing the conformal coordinate ( = e k ^/k one writes the general solution of the bulk 
equations (15) as a combination of Bessel functions, 



with Ai, Bi, A 2 , B 2 being complex numbers. For each value of the mass m there are two 
eigenvectors Xm, r — 1,2, of Eq. (17), which are normalized as follows, 



The normalization factor is determined by the asymptotics of the modes at large (. Making 
use of the asymptotics of the expressions (19) one translates (20) into the relation between 
the coefficients A, B, 




Xi = AMmC) + B 1 N 2 (mC) , 
X2 = m( [A 2 Ji(mQ + B 2 Ni(m< 



0] 



(19a) 
(19b) 




(20) 



2k 



■rs 



(21) 



m 
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Let us turn to the boundary equations (16). In terms of the coefficients A, B they take the 
form, 

*(Mt)-s*(t)) + *( w '(t)-^(t)) 

-^'>(t)-*iMtH- ™ 
-<*(t)-*£M7) + Mt*(tK*(t)) 

+ *(tMtKMt)H- < 22b > 

It is convenient to choose the first eigenvector of the operator (17) to satisfy the relation 

AxS 1) +/ix? ) L =n = o. 



I 2 = 

Then Eqs. (22) yield 

A? = CO'iV, (12) , B<» = -C«7, (£) , 

41) = _c.(u « jv (?) , b<" = cm^j m , 

/im \ k J fim \ k J 

where the normalization constant is fixed by Eq. (21). At m <C k it has the form, 

implying the following amplitude of the graviton field on the brane 



Xi 



(i) 



" \2k 



These modes are analogous to the continuum graviton spectrum in the RS case. They are 
completely delocalized and provide corrections to the four-dimensional Einstein gravity only 
at short distances r ~ 1/k. 

The orthogonal mode is fixed by Eq. (21). In the general case the corresponding ex- 
pressions are rather complicated. A considerable simplification occurs, if one assumes the 
hierarchy between the parameters, A^C/i^Cz/^cA;~ M5, where M5 = (IQnG^ 1 ^ is 
the five-dimensional Planck mass. Since the physics at large distances is governed by light 
modes we are interested in the region of masses m 2 <C Xk. One obtains, 

n fi ( m \ 3 / 2 1 



1 v \2kJ \ 2k J ' 1 v \2k 



v \2kJ 1 _ 2A ln^ ' 

v k 



2 " 2 V m J 2 " 2v\m) 1 - ^ In 
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where C is the Euler constant. An important feature of this mode is the divergence of the 
amplitudes of the fields on the brane in the small mass limit, 

,f>l = .a(^y ! 

*' U -° v\m 1- 2* In? ' 

x ' V k 

( 2)| A (2kV' 2 1 

x ' v k 

(2) 

Thus, the modes give dominant contribution to the Green's function of the operator 
(17) at large distances on the brane. Their collection comprises the quasilocalized four- 
dimensional graviton. 

To make the last statement more precise, let us consider production of gravitational waves 
by an external periodic transverse traceless source on the brane, 

T%*{x,z) = 6(z)e-**f ^ 3 e ipx d%\p)T^(p). 

The resulting gravitational field is given by the convolution of the source with the Green's 
function 

G(x - x'; lu) = -IQitG 5 J dr G u (x, x'; z = z = 0)e~ iuJT , (23) 

where r = t — t', and G n ( ) is the upper left element of the retarded Green's function 

of the operator (17), 



roc A 



(24) 



A/-1 POO 2 

G(x - x'; uj) = *- / dm V ( xS(0) ) e*" , (25) 



Substituting Eq. (24) into (23) one obtains, 

roc 

/ rfm 



where r = |x — x'|, p w = \J uj 2 — m 2 when m < uj and p w = i\Jm 2 — to 2 when m > a;. 
We see that the gravitational field on the brane has the form of a superposition of massive 
four- dimensional modes. Only modes with m < uj are actually radiated, the other ones 
exponentially fall off from the source. Thus, as long as we are interested in the gravitational 
waves, we can integrate in (25) only up to m = uj. This means, in particular, that at uj <C k 
one can neglect the contribution of the modes Xm m Eq. (25). Let us introduce 



m, 



= ke~ u ' 2X , (26) 
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and study the regime m c <C u <C k. The Green's function (25) takes the form 

--/o -(1-^lnf) 2 ' ( } 



where = G 5 k is the four-dimensional Newton constant. At r u/m 2 c the integral in 
Eq. (27) is saturated by m ~ m c and we obtain the usual four- dimensional expression for 
the gravitational wave. We stress that this wave is a superposition of massive modes with 
masses of order m c . In the opposite limit r uj/m 2 c the integral is damped by the rapidly 
oscillating exponent, so one obtains 

nt i \ ^Gn v ■ 

G(x - x ; a;) = nr- e . 

r Aln^ 

The gravitational wave gradually dissipates into the fifth dimension. 

A few comments are in order. First, the graviton mass scale m c turns out to be indepen- 
dent of the parameters Xi, x 2 in the action (1): 

v 1 
2A ~ l6irG 5 v 2 ' 

This is not very surprising. Indeed, from the very beginning we have been interested in the 
modes with masses much smaller than the energy scale associated with x x and x 2 . Then, 
the potential term for vector fields is effectively frozen at zero and its parameters do not 
affect dynamics. Rather, the condition that the potential must be zero, generates mixing 
between gravitons and vectors which leads to quasilocalization of graviton. The strength of 
mixing at low energies depends only on v, and so does the graviton mass. 

Second, very small values of m c are generated without fine-tuning. For example, taking 
k, M 5 as large as the Planck mass and v w (M5/5) 3 / 2 one obtains m c as small as (10 28 cm) _1 , 
that is the inverse of the present horizon size of the Universe. 

Finally, though the Lorentz symmetry is broken in the background (2), the Green's 
function (24) for the tensor modes has Lorentz-invariant form. 1 The place where the Lorentz- 
breaking is essential is that the very notion of the tensor modes refers to the preferred 
reference frame. 



1 This allows to interpret the additional enhancement of the distance r <~ uo/m? c , at which the dissipation 
of the gravitational waves sets in, by the factor uj/m c as compared to to" 1 , as a relativistic effect, see Ref. [2]. 
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4 Vector modes 



The gauge conditions (3), (6) result in the following equations for the vector perturbations, 

po^i + P 2 ^2 = , (28a) 
Po^s + P 2 ^4 + P 2 ^5 = , (28b) 

where we again omitted the polarization index a. From Eqs. (4), (5) one obtains the bulk 
equations for the vector modes, 

^-2^7 + 1^ = 0, 7 = 1,2, (29a) 
v 2 

-< + 2^-^/ = 0, 7 = 3,4,5, z>0, (29b) 

The junction conditions on the brane read: 

V>i + 2A;^i = , (30a) 

^ 2 + 2hh = \(ih + — + —) , (30b) 
\ v v J 

4>' 3 = ^ = , (30c) 

^. xV L + ^ + *) . (30d) 
\ v v J 

Equations (30a), (30b) together with the gauge fixing condition (28a) imply 

2^4 



^2 + 



+ ^ = 0. (31) 

z=0 V 



Thus, the system (30) reduces to a set of homogeneous boundary conditions for the vector 
part of the gravitational perturbations and for the free vectors in the RS background. In 
other words, mixing between gravitational and vector fields has no effect in the sector of 
vector modes. The relations (28), (31) imply three independent vector modes, parametrized, 
say, by the functions ^2(2), ^3 (z), ip4,(z). 

Let us show the absence of physical zero modes in the vector sector. Assume that there 
is a normalizable mode with p 2 = 0. Then, from Eqs. (29) we find that 1^2 = C2e _2fc ' z ', 
and that ^3, ^4 do not depend on z. The latter fact means that the corresponding vector 
a£ is independent of z, see Eq. (14e). Consequently the vector field strength 7 1 " is also 
independent of z, and, for the mode to be normalizable, F£ v must be zero. This implies that 
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the vector field of the zero mode is pure gauge in the bulk. From Eqs. (14b), (14c) we 
find that the metric of the zero mode is also pure gauge, 

h, u (x, z) = + U»>~ 2m , ^ = -iC^e^e-** . 

Thus, in the bulk, this mode can be removed by residual gauge transformation left after 
fixing the transverse-traceless gauge (3), (6). The only degree of freedom we are left with is 
the vector component vp e of the field f3 a , see Eq. (14f). This field resides on the brane and 
describes explicit breaking of the [?7(1)] 3 gauge symmetry by the brane action, cf. Eq. (12). 
But ip 6 vanishes due to Eq. (31). Hence, all components of the mode vanish: there is no 
normalizable zero mode at all. One concludes that there is only the continuum spectrum 
in the vector sector with p 2 > 0. The corresponding modes are completely delocalized and 
have the same form as in the pure RS case. 

It is worth mentioning that at p 2 = the mode ipQ vanishes due to field equations and 
boundary conditions. This could be dangerous, if this mode became propagating in a non- 
trivial background. As we show in Appendix B, this is not the case, which implies also that 
ipQ remains non-dynamical to higher orders in classical perturbation theory. 



5 Scalar sector 

We now turn to the scalar sector. Having in mind the issue of vDVZ discontinuity, we will 
consider the gravitational field produced on the brane by an external energy-momentum 
tensor T^ 1 localized on the brane. We do not restrict our analysis to static sources and 
consider external energy-momentum tensor depending both on time and space. Passing to 
the 4-dimensional Fourier representation we decompose it as follows, 

7S* = ti, T™ t =p i t 2 + ... , TV* = ^5 ij -^jt 3 +p i p j t 4 + ... , 

where we wrote down explicitly only the scalar components, which are of interest to us. The 
energy-momentum conservation implies 

Poh + p 2 t 2 = o , p t 2 + P % = o . 

Below we will consider two gauge-invariant scalar potentials, 03 and <§ = 0i + 2p o 02 +Po^4> 
of the metric produced by the external source on the brane. The quantities with the bar 
are defined via decomposition of the metric on the brane in accordance with (14a)-(14c). 
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The potentials 03, $ are related to the standard Bardeen variables $, \I/ (in the notations of 
Ref. [26]) as follows, 3 = 2\1>, $ = 2$ + Using the relation (11) between the metric 

induced on the brane and the bulk metric we obtain, 

3 = 3 | 2=o + 2fe, (32) 

2kv 2 

$ = (0! + 2p o 02 + Po04) + — r e . (33) 

z=0 p 2 

Our purpose is to compute the potentials 3 and $ and to compare them to the expressions 
obtained in the usual four-dimensional Einstein gravity. 

In the bulk, the fields obey the following equations, see Eqs. (4), (5), 

- 2k 2 ^ + = , 1 = 1,2,3,4, (34a) 

2 

- </>'! + 2k<p'r - ^0/ = 0, 1 = 5,6,7, z>0 , (34b) 

The gauge conditions (3), (6) take the form, 

0x - 203 - p 2 04 = , (35a) 

Po0i + P 2 02 = , (35b) 

PO02 + P 2 04 = , (35c) 

PO05 + P 2 07 = . (35d) 

As described in Appendix A, the junction conditions in the presence of the external energy- 
momentum tensor on the brane can be cast into the following form, 



<p[ + 2k(j)i = 2p 2 e + 87rG 5 ti , 


(36a) 


02 + 2k(p 2 = -2p e + 8nG 5 t 2 , 


(36b) 


2 

03 + 2/c03 = -p 2 e - ArrG^h , 

P 2 


(36c) 


2v 2 

<f)' 4 + 2k<f) 4 = -^-e + 8ttG 5 U , 
p 2 


(36d) 


05 = 0, 


(36e) 


06 = -T^vph - 2nG 5 ^v{-t 1 + 2t 3 + p 2 t 4 ) , 


(36f) 


0' 7 = , 


(36g) 


- p 2 e = X lP I 03 + -0 6 + 2ke + —^(-ti + 2t 3 + PU) , 


(36h) 
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where 



2„,2 



V\ = 2x 1 i> 



P = 



Ai + A 2 



3Ai + 2A 2 

These equations should be supplemented by Eq. (51a) in Appendix A, which determines the 
longitudinal component 8 in terms of the other fields. Note that Eqs. (36c), (36f), (36h) 
form a closed system. Once the solution of this subsystem is found, the other equations are 
solved in a straightforward manner. 

Before proceeding with the calculation of the gravitational field produced by the source, 
let us make sure that the scalar sector does not contain instabilities. So, one temporarily 
sets the external source equal to zero in Eqs. (36). First, we show that there is no tachyonic 
mode, which would correspond to negative or complex p 2 . For such a mode, let us introduce 
w, such that k 2 w 2 = —p 2 and Rew > 0. Then, for a normalizable mode, one would have 

03 = U 3 K 2 (we k ^) , 6 = UeweWK^weW) , 

where K ± , K 2 are modified Bessel functions. From Eqs. (36c), (36f) we obtain 

kwe rr 3vivke 



Ua 



Inserting these expressions into Eq. (36h) one obtains the following relation, 

Aip 



w + 



k 



K (w) 3i/i Ki(w) 



. 



(37) 



Ki(w) Ai K (w)_ 

Let us show that the real part of the expression in square brackets is positive in the right 
half-plane Rew > 0. Indeed, on the boundary of this half-plane one has, 



Re 



Re 



K (w) 

KiH 
K (w) 



w=iy 



n\y\(J 2 (y) + N 2 (y)) 
= 1 > . 



>0, 



I »oo,— 7r/2<Arg«)< 7r /2 



Hence, Re 



K (w) 



, being a harmonic function, is positive everywhere inside the half-plane 



K (w) 



>0, 



Rew > 0. This ensures the positivity of the real part of the inverse function, Re 
and thus of the whole expression in square brackets in Eq. (37). One concludes that Eq. (37) 
has no solutions in the right half-plane, implying the absence of tachyonic modes in the 
scalar sector. 

Second, let us demonstrate that the model is also free of ghosts. The ghost mode, if 
any, must be localized on the brane. Indeed, the normalization of the modes of contin- 
uum spectrum is determined entirely by the bulk action, which by itself is free of ghosts. 
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It is straightforward to see that the modes corresponding to strictly positive p 2 belong to 
continuum part of the spectrum. The only dangerous eigenvalue is p 2 = 0. But the corre- 
sponding mode is unphysical. Indeed, for this mode 02 = t/2e _2fc ' z ' and from Eq (36b) one 
obtains e — 0. Substituting p 2 = p^ into the gauge fixing conditions (35b), (35c) one finds 
X = _ Po [/ 2 e- 2fe W, 4 = _g e - 2fe l< Then, Eq. (35a) yields 3 = 0. It follows from Eq. (34b) 
that the scalar component 07 of the vector fields does not depend on z. From Eq. (35d) we 
find 05 = — po07- Finally, from Eq. (36h) we obtain 6 = 0. The surviving mode is pure 
gauge in the bulk, see Eqs. (14), 

The bulk fields can be removed by a residual gauge transformation, leaving the scalar com- 
ponent 08 (see Eq. (14f)) of the fields f3 a on the brane. But the latter is zero according to 
the field equation (51a). Thus, all the components of the zero mode vanish. There is only 
continuum spectrum of modes with p 2 > in the scalar sector. 

We show in Appendix B that the mode 08 remains non-dynamical in a non-trivial back- 
ground as well. Hence, no new propagating degrees of freedom appear both in non-trivial 
backgrounds and in higher orders of classical perturbation theory. 

Having established the absence of instabilities, we now return to the field produced by 
the source. Our strategy is to express the functions 03 (2), 4>&{z) in terms of the external 
source and the brane bending e using the bulk equations with the boundary conditions (36c), 
(36f), and to insert the result into Eq. (36h). This is most easily done by combining the bulk 
equations with the boundary conditions into the Schrodinger type operators and introducing 
the Green's functions of these operators. The latter satisfy the following equations, 

-u 2 G$(z;p) + 4k 2 u 2 G 3 (z;p) - Ak5(z)G 3 (z;p) -p 2 G 3 (z;p) = S(z) , 
-u 2 Gl{z- p) + 2ku 2 G' 6 (z; p) - p 2 G 6 (z; p) = 8{z) . 

Imposing the radiation (outgoing wave) boundary conditions at z — > ±00 one obtains, 

G3{ZtP) = ~2pH[%/k) ' (38) 

G 6 {z,p) = — , (39) 

2pH^(p/k) 
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where p = a/f> Rep > 0. Making use of these Green's functions we find, 



.(0) 



H?{p/k) 
pH^\p/k) V2Ai 



e + 47rG 5 



ti 



3ui 



vp 2 e + 2nG 5 —v(-t 1 + 2t 3 + p%) 



(40) 
(41) 



Substitution of these expressions into Eq. (36h) yields, 



pH$\p/k) 



P 2 + Xip 



H?{p/k) 



3z/ip 



pH[ 1 \p/k) 



H$\p/k) 



4ttG, 



-T 



H?{p/k) 
pH^ip/k) 



3z/i/> 



pH^ip/k) h 



where we have used the relation 

2k 



pH^ip/k) _pH^{p/k) 



47rG 5 \ lP 



(42) 



H?{p/k) P 2 



= 



H^ip/k) " H^ip/k) 

and introduced the notation T = t\ — 2t 3 — p 2 t 4 for the trace of the energy-momentum tensor. 

Further analysis depends on the value of the parameter x 2 . Let us first consider the case 
x\ = 0, which in Eq. (1) corresponds to the absence of the first term in the potential for 



vector fields. 2 In this case Ai = v\ 



0, and Eq. (42) reduces to, 
AnG 5 , 



e = 



3p 2 



-T . 



(43) 



The last expression coincides with the result obtained in the pure RS case, cf. [24]. Substi- 
tution of (43) into Eqs. (40), (41) leads 3 to vanishing of the scalar part of the vector fields, 
6 = 0, while the scalar components of the metric are the same as in the pure RS case. 
Hence, there are only short- distance corrections to the scalar part of the metric in the case 
k\ = 0. At distances r 3> 1/fc, i.e. at small momenta p <C k, one arrives at the same 
expressions as in the four- dimensional Einstein gravity, 

87tGn 



-SttGW 16tiG n + 

(P 2 ) 2 P 2 



(44) 



where Gn = G$k is the four-dimensional Newton constant. Evidently, the gravitational field 

is free from the vDVZ discontinuity. In particular, the Newton law — the field of a static 

source — does not get modified at large distances in the case h\ = 0. 

2 Wc do not know whether vanishing of the parameter k{ can be ensured by any symmetry requirement. 
3 Note that though Ai = v\ = 0, the ratio ^ = 8jv q 5V -i in (41) is finite. 



18 



Another possibility is that the parameter x\ is large. In this case we are interested in 
the behavior of the fields at distances r ^> 1/Vi, i.e., at small values of the momentum, 
p <C fc, v\. We also assume as in Sec. 3 the hierarchy Ai Ci/i. The dominant contributions 
to the expressions in the square brackets in Eq. (42) come from the terms proportional to 
H^~\p / k) / Hq~\p / k) . It is worth noting that these terms come from the scalar component 
06 of the vector fields in Eq. (36h). Keeping only the leading contributions at low momenta, 
one obtains from Eq. (42), 

£ = ^ T + ^ tl V 1 ln k- (45) 

This expression for the brane bending differs from the formula (43) in the pure RS case 
by the last term, which is small at moderate momenta. However, it becomes important 
and is logarithmically large relative to the first term in the far infrared. Making use of the 
expression (45) for the brane bending e and the Green's function (38) one determines the 
rest of the gravitational modes, 0i, 02, 04- Finally, these expressions, being inserted into 
Eqs. (32), (33), yield the scalar components of the metric induced on the brane, 

03 = -SttGaA , (46) 
P 

- 8nG N p 2 16ttG n \QnG N p 2 X 1 p 

(p 2 ) 2 p 2 (p ) vx k 

The expression (46) coincides with the result of the usual four-dimensional Einstein gravity. 
Thus, the potential 03 does not get modified at large distances at all. As to the second 
potential, the first two terms in Eq. (47) also coincide with the result of conventional gravity. 
The last term provides the long-distance correction, which becomes important at the distance 

~~ k ~ k 

This distance coincides with l/m c , where m c is the graviton mass scale (26). 

To make the physical consequences of the formula (47) more clear, let us calculate the 
gravitational field produced by a point-like static source. We take T 0Q = M5(x), T 0i = = 
0. Setting t\ = M, t 3 = 0, p 2 = — p 2 in Eq. (47) and performing Fourier transform, one 
obtains, 

= , x 2G N M / 2Ai , , \ 

$(r) = — 1 -Infer . (48) 

r \ V! J 

while for the other gauge invariant potential we have, 

M r ) = • ( 49 ) 
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These expressions are valid up to corrections at small distances r ~ \ jk,\jv\. The potential 
(48) is the analog of the Newton potential in our model, it is responsible for gravitational 
interaction between nonrelativistic massive objects. Remarkably, a contribution to $(r) ap- 
pears, which grows logarithmically with the distance as compared to the standard expression 
{—IG^Mjr). Even more strikingly, the sign of this contribution is opposite to that of the 
standard expression. As a result, the expression (48) describes gradual weakening of gravi- 
tational attraction with the distance; at large distances, r > l/m c , attraction gets replaced 
by repulsion. Thus, our model provides an example of a ghost-free theory with antigravity 
at ultra-large distances. 

At first sight it seems surprising that in the model with quasilocalized gravitons we 
obtain a contribution to the Newton potential which grows at large distances, as compared 
to the standard four-dimensional expression. Following the conventional line of reasoning 
one could infer that, as gravitons dissipate at large distances into the fifth dimension, the 
gravitational interaction should become weaker at large distances, as compared to the four- 
dimensional case. However, this line of reasoning is incorrect. In theories with Lorentz 
symmetry breaking, the potential produced by an external source is not directly related to 
the spectrum of propagating degrees of freedom. This point is illustrated by four- dimensional 
Lorentz violating massive electrodynamics, considered in Refs. [27]. In that case the electric 
potential of static sources falls off as 1/r in spite of the fact that all propagating modes are 
massive. 

The origin of the logarithmically enhanced antigravity in our model can be understood 
as follows. A point mass gives rise to perturbations of the vector fields which interact with 
matter via mixing with the metric. The gravitational field is produced by the total energy- 
momentum tensor composed of Tf^ and the energy-momentum tensor TY V of the vectors, 
given by Eq. (9). The latter tensor falls off slowly from the localized external source, and 
dominates at large distances. Our results indicate that, insofar as the Newton potential is 
concerned, the vector fields mimic the effect of negative energy. This is possible because T^ v 
violates the weak energy condition. 

The gravitational potentials (46), (47) are free from the vDVZ discontinuity. Indeed, in 
the limit of vanishing graviton mass scale m c , which corresponds to \\jv\ — > 0, one recovers 
the usual four dimensional expressions (44). However, from the phenomenological point of 
view the difference between the two gravitational potentials (48) and (49) results in a severe 
phenomenological constraint on the value of the parameter \\jv\. Measurements of the light 
deflection by the gravitational field of the Sun [28] require \\jv\ < 1CT 5 . This value is 
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not unnaturally small: it corresponds to v ~ (0.027M 5 ) 3/2 . However, it pushes the range 
l/m c , where the antigravity sets in, far beyond the present horizon size of the Universe. 
The mass scale m c in this case is very small and it is not clear whether it can significantly 
affect the physics within the present horizon. We discuss possible ways of avoiding this 
phenomenological constraint in the concluding section. 

6 Discussion 

In this paper we presented a Lorentz-violating brane-world model, where gravitons are not 
completely localized, but are rather quasilocalized on the brane. In other words, the four- 
dimensional graviton is a collection of KK modes from the continuum spectrum. The char- 
acteristic mass m c of these modes is exponentially small when expressed in terms of the 
parameters of the Lagrangian. We demonstrated that the model is free from tachyonic 
and ghost-like instabilities. We calculated the metrics produced by an external energy- 
momentum tensor on the brane, and found that the model does not suffer from the vDVZ 
discontinuity. 

The key observation behind the model is that while the gravitational perturbations in 
the Randall-Sundrum setup contain a zero mode localized on the brane, the bulk vector 
fields are completely delocalized. Thus, mixing between the vectors and the would-be zero 
gravitational mode forces the latter to dissipate into the fifth dimension. This mixing between 
tensors and vectors becomes possible when the Lorentz symmetry is spontaneously broken 
by non-zero VEVs of the vector fields. These considerations are rather generic. We believe 
them to be applicable to a broad class of generalizations of the Randal-Sundrum model with 
spontaneous Lorentz symmetry breaking by bulk fields. In particular, the vector fields can 
be replaced by form-fields of higher degrees, which are also completely delocalized in the 
Randall-Sundrum background. 

For general form of the potential term for the vector fields on the brane the model ex- 
hibits antigravity at ultra-large distances. Namely, the structure of the Newton potential is 
such that gravitational attraction between nonrelativistic massive objects gradually weakens 
as the distance increases, and gets replaced by repulsion at r > l/m c . While this property 
is theoretically appealing, it puts severe phenomenological constraints on the model. The 
reason is that the second Bardeen potential, characterizing spatial part of the metric of a 
static source, does not get modified at large distances and thus differs from the Newton po- 
tential. Observations of light deflection by the Sun constrain the relative difference between 
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the two gravitational potentials at the level of 10~ 4 . This translates into the constraint 
m c < Mpi exp (— 10 5 ). This value is negligible compared to the present Hubble parameter of 
the Universe. It is doubtful whether the model with so tiny graviton mass scale can lead to 
any interesting phenomenology, in particular, to accelerating expansion of the Universe at 
the present epoch. 

One possibility to avoid this phenomenological problem is to fine-tune the parameter K\ 
of the vector potential in the action (1) to zero. Then, both Bardeen potentials have at large 
distances the same form as in the four-dimensional linearized Einstein gravity, while graviton 
is still quasilocalized and gravity waves dissipate into extra dimensions. The graviton mass 
scale m c in this case can be comparable to the Hubble parameter or even larger. A drawback 
of this approach is fine-tuning of X\. We are not aware whether vanishing of this parameter 
can be imposed by any symmetry. 

Another way out is a generalization of the model considered in this paper. Let us sketch 
a particular example. One introduces a scalar field with dilaton-like coupling to the vector 
fields on the brane. Namely, one replaces the combination g^ v A a ^A b v in the brane action in 
Eq. (1) by g^ v e aip A a ^A b u . The modified model allows for the Lorentz- violating background 
(2) with tp — 0. The dilaton does not affect tensor and vector sectors of the linearized 
perturbation above this background, so the modified model also incorporates quasilocalized 
gravitons. On the other hand, relative difference between Bardeen potentials depends on the 
dilaton coupling, and is essentially proportional to 1/a 2 when a is large. Thus, the constraint 
from light deflection is satisfied once a > 100. We will report more on phenomenology of 
our model and its generalizations elsewhere. 

Two other important issues which we leave for future investigations are the quantum 
structure of the theory, and cosmology in the model (1) and its generalizations. Of special 
interest is at what scale the strong coupling at the quantum level sets in, and what kind of 
late-time cosmological evolution is obtained in these models, in particular, whether they can 
account for the cosmic acceleration at the present epoch. 
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A Junction conditions in the scalar sector 

In this Appendix we consider the junction conditions for the scalar modes in the presence of 
an external energy-momentum tensor on the brane. The external source should be added to 
the r.h.s. of Eqs. (13) according to Eq. (7). As a result one obtains the following junction 
conditions on the brane, 

A / 4 2d 2 d 2 2in 2 

(/)[ + 2A;0! - 2p 2 e =-± 20 3 + p 2 4 + -06 + —07 - K~£ + 6fce + 



6 \ V V K V 

+ 87rG 5 (\t x + \t z + , (50a) 



S 3 ° 3 

02 + 2A;0 2 + 2p e =8nG 5 t 2 , (50b) 



/l n 2 2 2n 2 n 2 2?'n 2 \ 

4 + 2*. = A(-0 3 - ^ 4 + ^« - + |, - J» , (50c) 

/ 2 2 4 4 2 4i 

0^ + 2/c0 4 - 2e = A I - ^03 + -04 - 7^">6 + —07 ~ 77^ + 7^08 

+ 8 ^ 5 (^-^ 3 + ^)' (50d) 



05 =0, 



(50e) 



? ■' 2 ' />/ - 4 - . J/, ) ^ '(^03 -^04-7?" ^ P 



= f I (/> 3 ~| c/^6 



^4 T — V6 

3 3v 



3f Vl ~^ 3k 



2ip 2 
~3?T 



8 , 



0' 7 =X 2 W 3 f 4 + 



2 , 1 2fc 2i 
-07 - T e + + —< 
v k p z v 



, ... 2 , 2 4 
4 W| -V 03+ 3 04 -3p^ ( 



4 J 2 4i 
»6 4- t^-07 - —£ 4 



3w 



3A; 3v 



(50f) 



(50g) 
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Combining these equations with the gauge fixing conditions (35), we obtain 



Ai -P 2 04- 



+ A 5 



2p 2 



P 2 rt , 2ip 2 
+ Z-e - Ike — « 

k v 



2=0 

2 , 2p 2 , 4 

^03 - ^-04 + — < 

3 3 3v 



Ap 2 
~3v 



2p 2 4ip 



2p 2 £ =^ <! ( 203 + P"04 + -06 4 



2 P 2 



2=0 

2 



3A; 



3w 



2 = 



p- 2*p 2 
— £ + 6ke H S 

K V 



(51a) 



(51b) 



The first of these equations plays the role similar to that of Eq. (31) in the vector sector: 
it determines the field 8 in terms of the other fields. Substitution of Eqs. (51) back into 
Eqs. (50) yields the system (36) considered in the main text. 



B Absence of ghosts above a non-trivial background 

We saw in the main text that the modes 0s and -06a, accounting for the explicit [?7(1)] 3 
gauge symmetry breaking on the brane, vanish the the linear order in perturbation theory 
above the background (2). These vanishing modes are generically dangerous. In higher 
orders of classical perturbation theory (i.e. above non-trivial backgrounds) kinetic terms for 
these modes may arise, rendering them propagating. Depending on the sign of kinetic term 
some of these modes may become ghosts. To study this issue we investigate the ultraviolet 
behavior of the perturbations of the vector fields above a non-trivial vector background 
close to the static background Af = vSf. For the sake of simplicity we neglect gravity 
perturbations. This can be done consistently by taking the limit G 5 0; we believe that 
switching on gravity does not spoil our results. The non-trivial background is chosen to be 
locally space-like = (0, A a ). In this Appendix we show that the equations of motion 
for the linear perturbations above the non-trivial vector background imply vanishing of the 
dangerous modes, in complete analogy to the situation above the static background (2). 
This demonstrates the ultraviolet stability of the model (1) to higher orders in classical 
perturbation theory. 

Let us consider the brane part of the quadratic action for the fluctuations a a of the vector 
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fields above the non-trivial vector background, 

S( 2 ) = J d 4 x -2 (x 2 + x 2 ) [A? b A c p a v b a c u + A» b a c (i a v b A c u + A» b a c ^a u c A b u ) 

2x 2 i 
- 2x x V ■ a u b a b v - — ± [A? b A\a v c < + 2A» b a\A i c a c ^ , 

where summation over internal indices b, c = 1,...,3 is assumed, and the metric on the 
brane is taken to be flat. The bulk part of the action is irrelevant for us as we are interested 
in longitudinal modes = d^f3 a which are pure gauge in the bulk. 

Let us consider equations for the spatial components of the fluctuations, a" = dif3 a . It 
is convenient to represent them in the matrix form introducing 3x3 matrices a = a® and 
A — Af. By a suitable transformation, the background matrix A can be made symmetric, 
A = A T . Then the equation of motion for the fluctuations reads, 

4 (x 2 + x 2 ) [A 2 d + Aa T A + ai 2 ) + 4xjV • a + ^ (^2Tr[ai] • A + Tr[i 2 ] • a) = . (52) 
The symmetric part of the fluctuations, s" = dif3 a + d a (3\ obeys the homogeneous equation, 

A[s] = , 

where the linear operator A can be read from Eq. (52). For the static background (2b) the 
operator A is non-degenerate. Thus, it is non-degenerate for non-trivial backgrounds, which 
are close enough to (2b). One concludes that the matrix s vanishes, 

di(3 b + dbft = . (53) 

Differentiating Eq. (53) twice yields the Laplace equation for the divergence, Ad a f3 a = 0. 
Imposing vanishing boundary conditions at spatial infinity we obtain d a f3 a = 0. Then, 
differentiating Eq. (53) once, we find that all the components f3 a vanish. 

Thus the longitudinal modes 0g and ipQ a do not become propagating in non-trivial back- 
ground. This result implies the absence of ghosts and/or tachyons in the model (1) to higher 
orders in classical perturbation theory. 
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